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Abstract 

We show that a special choice of light-cone gauge can greatly sim- 
plify the calculation of the classical color field created in the initial stages 
of nucleus-nucleus collisions. Within this gauge, we can in particular 
construct explicitly the conserved color current and calculate exactly the 
gauge field immediately after the collision. This field is used as a bound- 
ary condition in an iterative solution of the Yang-Mills equations in the 
forward light-cone. In leading order, which corresponds to a linearization 
of the Yang-Mills equation in the forward light-cone, we obtain a simple 
formula for the spectrum of gluons produced in nucleus-nucleus collisions. 
This formula reproduces exactly the known formula for proton-nucleus 
collisions, where kt— factorization is recovered, while the latter property 
apparently breaks down in the case of nucleus-nucleus collisions. 
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1 Introduction 



It is believed that strong color fields are created in the initial stages of nucleus- 
nucleus collisions at high energy. In the color glass formalism, or its early 
formulation by McLerran and Venugopalan [1] these fields are determined by 
solving the Yang-Mills equations for given distributions of color charges carried 
by the two nuclei. An average over these color charges is performed in the 
calculation of the observables. In the original formulation, the distribution of 
the color charges carried by the nuclei was assumed to be Gaussian. However, 
the definition of the color charges involves a separation of scales between the 
short wavelength partons, treated as sources, and the long wavelength ones, 
treated as fields. Independence of the final results upon this dividing scale leads 
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to a renormalization group equation describing the non linear evolution of the 
parton density in various forms [2-9]. 

We shall not be concerned in this paper with this non linear evolution (for 
reviews see e.g. [10,11]), but shall focus merely on the determination of the 
initial classical field produced by a given distribution of color charges. This 
problem has so far resisted a complete analytical treatment. Various numerical 
studies have been performed [12-14], based in particular on Refs. [15,16] (sec 
also [17]), and several approaches have been suggested: Kovchegov's treatment 
relies on the conjecture that final state interactions do not play any role, in ad- 
dition to other assumptions [18]. More recently, Balitsky proposed a symmetric 
expansion in powers of commutators of Wilson lines [19]. The only case for 
which a complete analytic treatment exists is that where one of the two sources 
can be considered as a weak color source. This is the situation for instance in 
proton-nucleus collisions. One can then linearize the Yang-Mills equations with 
respect the weak source of the proton, while treating all the high density effects 
in the strong source of the nucleus [20-25] . This allows one to express the cross- 
section for gluon production as a convolution in transverse momentum space 
of the gluon distributions of the proton and the nucleus, a property known as 
k t — factorization . 

In this paper we report on progress in the analytic calculation of the color 
field produced in nucleus-nucleus collisions, and the resulting production of glu- 
ons. We show that the complexity of the calculation is greatly reduced by an 
appropriate choice of gauge. We choose an axial gauge that makes the field 
before the collision entirely concentrated in a narrow strip along the light cone. 
In particular, pure gauge fields that remain present in the forward light cone 
with most gauge choices, are absent in the present formulation. This work gen- 
eralizes to the case of nucleus-nucleus collision the study performed for the case 
of proton-nucleus collisions using a similar gauge choice [25] . 

The paper is organized as follows. In Section 2, we study the Yang-Mills 
equations in the light-cone gauge and show that important simplifications can be 
achieved by exploiting judiciously the remaining gauge freedom. We construct 
explicitly the conserved current and calculate exactly the field immediately after 
the collision. Then we propose an iterative solution of the Yang-Mills equations 
in the forward light-cone. The leading order involves a linearization of the 
equations and is used in Section 3 to calculate the spectrum of gluons produced 
in nucleus- nucleus collisions. A compact formula is obtained for this spectrum, 
which does not exhibit /cf-factorization. The latter property is recovered in the 
case of proton- nucleus collisions for which our formula yields the exact result. 
In Appendix C we show that our results can be recovered by using a variant of 
the Fock-Schwinger gauge. 
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2 Solving the Yang-Mills equations in light-cone 
gauge 

We consider a collision of two nuclei, denoted A and B, that are moving at 
nearly the speed of light, respectively in the —z and the +z directions, and 
carrying color charge densities p A and p B . Because of the Lorentz contraction, 
these densities have a small longitudinal extent, which translates into small 
extents, of order e, in the x + direction for the nucleus A and in the x~ direction 
for the nucleus B, where x^ = (t ± z)/v2 denote the light cone variables (see 
Fig. 1). Thus p A and p B are functions of x + and x~ , respectively p A {x + ,x) 
and p B (x~ ,x), with x the transverse coordinates. 

Note that the integrals of p A over x + , and p B over x~ , are finite; it follows 
that the color charge densities diverge as 1/e when e — > 0. One often takes the 
limit e — > from the outset, writing p A B oc S(x^). However, the longitudinal 
structure of the sources plays a role in the average over the color sources [26,27] 
and also, as we shall see, in the analysis of the Yang-Mills equations, in particular 
in the determination of the continuity conditions of the fields at the boundary 
of the nuclei. We shall therefore keep e finite and take the high energy limit 
(e — > 0) only at the end of the calculations. 

The classical gauge field is solution of the Yang- Mills equations 

where = — igA^T c = — igA^ ■ T is the covariant derivative, T c is 
a generator of SU(3) in the adjoint representation, and A^, and are 
components of color vectors (of the adjoint representation). Equations such as 
Eq. (2.1) are to be read as equations between color vectors (whose color index is 
omitted to alleviate the notation). The color current satisfies the conservation 
equation 

D^ = 0. (2.2) 

In a gauge transformation, 

A" ■ T • T)fit _ J-Q^Qt, (2.3) 

W 

with £1 a matrix in the adjoint representation of SU(3), the color current trans- 
forms as 

j m -> nj^. (2.4) 

Before the collision, the only components of the color current are 

J + (x) = p B (x~,x), J-(x) = p A (x+,x). (2.5) 

Note that the precise definition of the charge densities, and the color currents, 
involves a choice of gauge: we shall assume throughout that the color densities 
are those of the covariant gauge (see Sect. 2.1). 
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In most of this paper, we shall work in the light-cone gauge A + = 0. The 
Yang-Mills equations read then [25] 

-d + (d^) - ig{A l ■ T) d + A l = J+, (2.6) 
D~d + A~ - D l F l - = J~, (2.7) 
d + F- 1 + D-8+A 1 -D j F jl =J\ (2.8) 

Note that the first equation. Eq. (2.6), does not contain any derivative with 
respect to x + (d~ — d/dx + ). Therefore, it can be seen as a constraint that 
relates the various field components at the same x + . 

2.1 The gauge field before the collision 

Before the collision, i.e., for t < 0, the gauge field is obtained by solving the 
Yang-Mills equations for each of the projectiles separately We shall denote by 
A^ (A 1 * ) the field created by the nucleus A (B) in the absence of the nucleus B 
(A). 

Consider first the nucleus A. The gauge field that it produces is not com- 
pletely fixed by the gauge condition A + = 0. One can exploit the residual 
gauge freedom to choose the gauge field at x~~ < to be either longitudinal or 
transverse. 

The longitudinal solution is obtained from Eq. (2.7), by assuming that the 
component A~ does not depend on x~ , i.e., d + A~ = 0, and that the transverse 
components vanish, A 1 = [25]. One gets then 

A-(x)=$ A (x+,x), A' A =A+ = 0, (2.9) 

where 

-d]_$ A (x + ,x)=p A (x + ,x). (2.10) 

With this gauge choice, the field created by the nucleus A before the collision 
has support in a small strip of width e along the semi-axis x~ < (see Fig. 1 
below). The corresponding current created by the nucleus A is simply 

J»{x)=5»-p A {x + ,x). (2.11) 

It is conserved since 

D^ A = (d+ - igA+ ■ T) J- = d+J- = 0, (2.12) 

where we have used successively the gauge condition = and the fact that 
J~ is independent of x~. 

The transverse solution can be obtained by performing a gauge transforma- 
tion on (2.9) that eliminates A~ . Denoting by U* the matrix that realizes the 
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gauge transformation, we have 

A + ■ T — > - — U f d + U = 0, 
ig 

A~ ■ T C/ f ($ . - T)U- —U^d-U = 0, 

z.g 

A*. ■ T — > -—UWU. (2.13) 

ig 

We shall verify shortly that the solution of these equations is provided by the 
following Wilson line in the adjoint representation: 



U(x + , x) = Vexp 



HI 



,x+ 

/ dz + $ A (z + ,x) -T 

J — OO 



(2.14) 



The symbol V in Eq. (2.14), and throughout, denotes path ordering (along the 
integration variable, here z + ). Note that the integral in the Wilson line (2.14) 
gets contribution only within the support of p A , i.e., for < z + < e. Therefore, 
U becomes independent of x + when x + > e, and U(x + < 0) = 1. 

With U given by Eq. (2.14), the second equation (2.13) yields A" = 0. The 
first equation (2.13) ensures that the gauge condition A + = is satisfied since 
U does not depend on x~ . The only non vanishing components of A are, as 
announced, the transverse components A 1 , which we may rewrite as follows 

r x+ 

A\{x) = - dy+uHy + ,x)d l <i> A (y + ,x), (2.15) 



where we have used the following identity (valid for any matrix U of the adjoint 
representation) : 

U ] T a U = U ab T b = T h ul a . (2.16) 

Note that since U = 1 for x + < 0, the field A l A is non zero only in the region 
x + > 0. Note also that A l A is independent of x + when x + > e. It follows that 
A 1 becomes, in the limit e — > 0, a discontinuous function of x + , at x + = 0. 

A 

The color current corresponding to this transverse solution is related to the 
current (2.11) by the same gauge transformation as in Eq. (2.13): 

J»(x) -> UU»(x) = 8^U\x+,x)p A (x+ 1 x). (2.17) 

This equation makes explicit the dependence of the color density, and the as- 
sociated color current, on the choice of the gauge. Note that the longitudinal 
solution (2.9) is common to both the light-cone gauge A + = and the covariant 
gauge [24]: indeed, when x~ < 0, d^A^ = d + A^ = when A^ A satisfies (2.9). 
This is the reason for the simple form (2.11) of the current for the longitudinal 
solution. In most of this paper, we shall use the longitudinal solution (2.9) as 
the field produced by the nucleus A. 
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Let us turn now to the nucleus B. A priori there exists for B the same gauge 
freedom as for A. However the condition A + = forces us here to choose the 
transverse solution (sec however Sect. 2.5). With p B (x~ ,x) the color charge of 
the nucleus B (in covariant gauge), and $ B (x~ , x) the solution of (cf. Eq. (2.10)) 



wc obtain 



Al ■ T 



-d\<S> B {x ,x) = p B (x ,x), 



-V ] {x-,x)d i V{x~,x) and A~ 



AZ 



0, 



where V is a Wilson line along x : 



V(x , x) = V exp 



ig / dz $ B (z , x) ■ T 

J —oo 



More explicitly, 



A i B {x) = -f dy-V\y-,x)d i $ B {y-,x) 

J — OO 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



Note that since <f> B (x~) has support in the strip < x~ < e, A 1 (x) vanishes 
for x~ < and is independent of x~ for x~ > e: as e — > a discontinuity builds 
up in A l B {x~) 

The current of the nucleus B is related to p B by (cf. Eq. (2.17)) 



J^(x) = 6^+V\x- 7 x)p B (x- 7 x). 
This current is conserved since 

D»J» = (d- ig A- ■ T) J+ = d-J+ = 0, 



(2.22) 



(2.23) 



where we have used successively A B = and the fact that J+ (x) is independent 
of x + . 



2.2 The conserved current 

We have verified in the previous subsection that the currents and , carried 
respectively by the nuclei A and B, are separately conserved before the interac- 
tion. In this subsection we shall obtain the expression of the conserved current 
for the interacting system. To this aim we note first that the conservation law 
D^J^ = 0, Eq. (2.2), constrains only one component of the current. We shall 
exploit the remaining freedom to force the current to satisfy the extra relations: 

D+J-=0, D-J+=0, J l = Q. (2.24) 

A current satisfying (2.24) obviously satisfies (2.2). Note also that the relations 
(2.24) are compatible with the conservation laws (2.11) and (2.22) satisfied by 



6 





Nucleus B 



Nucleus A 



Figure 1: The gauge field before the collision, in the light-cone gauge A + = 0. 
The nucleus A produces a static longitudinal field localized in a strip of width 
e near the light cone along the x~ axis. The field of the nucleus B is transverse 
and exists in the half plane x~ > 0. ft is independent of x^ for x~ > e. 



the current before the interaction. Thanks to the gauge condition A + = 0, the 
first equation in (2.24) reduces to 



0. 



(2.25) 



This implies that the current J is conserved along the x direction, and is 
therefore given by the initial current at x~ — > — oo, namely 



J-{x) = J-{x + ,x), 

with J~ the current of the nucleus A alone, given in Eq. (2.11) 
equation, D~ J + = 0, is solved by 



J+(x) =W(x) Jt(x~,x), 



with 



W(x) = Pexp 



i ft 



dz + A-(z + ,x~,x) ■ T 

-OO 



(2.26) 
The second 

(2.27) 
(2.28) 



The Wilson line involves the component A" (x) of the full solution of the Yang- 
Mills equations, which is a priori not known. However, the current in 
Eq. (2.27), proportional to p B (x~,x), has support only in the narrow strip 
< x~ < e along the x + axis. Therefore, in the limit e — > 0, only the value of 
A~ along the x + axis is needed to evaluate W . In Appendix A we show that 
this is given by the field $ A created by the nucleus A alone: 



lim A (x^ 



(2.29) 
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Thus, in the limit e — > 0, the Wilson line W of Eq. (2.28) reduces to the Wilson 
line U(x + ,x) given in Eq. (2.14). 

It follows that, in the light-cone gauge A + = 0, the conserved current in a 
nucleus-nucleus collision can be written thus 

J»(x) = 8^-p A (x + ,x) + 5 fl+ U{x + ,x)V' ! (x-,x)p B (x-,x), (2.30) 

The first Wilson line acting on ps brings the density from the covariant gauge to 
the light-cone gauge A + = 0. The second Wislon line, U(x + , x) summarizes the 
effect of the interaction. The formula (2.30) shows that, in the light-cone gauge, 
the conserved current can be constructed explicitly in terms of the sources p A 
and p B (defined in covariant gauge). In the interaction zone a color rotation 
takes place, but this is determined by the gauge fields created by the individual 
nuclei alone, and these are known analytically. In contrast, in covariant gauge, 
the determination of the conserved current requires the solution of the Yang- 
Mills equation in the interaction zone and this can be done only approximately 
[24]. Note that the interaction leaves the component J~ unaffected (as a con- 
sequence of (2.26)), while J+ acquires an x + dependence through the gauge 
rotation. There is therefore an apparent dissymmetry in the present description 
that has, however, no consequence in the calculation of the gluon production, 
as we shall verify explicitly in Sect. 3. Had we chosen for the field of the nucleus 
A the transverse solution (2.15), one would have obtained a more symmetrical 
current J M — > Wj^ = S^~W p A + 5^ + V^ p B , similar to that obtained in the 
Fock-Schwingcr gauge x + A" + x~A + = 0; this leads to the same spectrum of 
produced gluons (see Appendix C for a more complete discussion of the Fock- 
Schwinger gauge). 

2.3 The gauge field immediately after the collision 

We now turn to the determination of the gauge field immediately after the 
collision, that is, along the lines x + = e or x~ = e (indicated by a thick (red) 
line in Fig. 2), in the limit where e — ► 0. As we shall see, this field can be 
determined analytically. It will serve as a boundary value for the solution of the 
Yang-Mills equations in the forward light-cone > e, that will be considered 
in the next subsections. 

Before we start our analysis, it is useful to recall that the densities p A {x + , x) 
and p B (x~,x) diverge as 1/e when x^ = e — > 0. The same holds for the fields 
$> A (x + ,x) and & B (x~ ,x), solutions of Eqs. (2.10) and (2.18), respectively. The 
transverse components of the field remain finite, however, as can be seen for 
instance on Eq. (2.15) for A % A : the divergence ~ 1/e of $ A is eliminated by 
the y + integration ~ e. We shall now proceed to an analysis of the Yang-Mills 
equations in the region x^ < e, keeping only the dominant terms as e — » 0. 

By integrating the Yang-Mills equation (2.8) over x~ from to e, at fixed 





Nucleus B 





Nucleus A 



Figure 2: The field A^ in the forward light-cone is obtained by solving the Yang- 
Mills equations with initial conditions given on the border of a strip of width e 
along the light-cone (indicated by the thick (red) line) . The field in this strip is 
determined analytically as a function of the source p A and p B in subsection 2.3. 

x + and x, we get 

28-A l (e) - d l [A-{e) - $ J - ig(A~(e) ■ T) A l (e) 

-ig [ dx-{A- -T) d + A\x~) - [ dx~ D j F jl ( X -) = 0. 
Jo Jo 

(2.31) 

Simplifications occur in the limit e — * 0: the second term drops out because 
of Eq. (2.29); the last term vanishes since Fij and its derivatives with respect 
to the transverse coordinates remain finite as e — ► 0; finally, using again the 
continuity of A~(x~), Eq. (2.29), one can perform the integral in the fourth 
term and show that its contribution is identical to that of the third one. Thus, 
Eq. (2.31) reduces to 

d-A t (x+,e,x) - ig[$ A (x + ,x) ■ T] A i (x + ,e,x) = 0. (2.32) 

This equation governs the x + dependence of the transverse components of the 
field, along the axis x~ = e. It is valid in fact for all x~ < e. It is easily solved, 
to yield: 

A l (x) = U{x + ,x)A i B (x-,x), x~<e, (2.33) 

with U given by Eq. (2.14) and A l B by Eq. (2.21). We have used the fact that 
A l {x+ = 0,x-,x) = A B (x~,x). 
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Similarly, by integrating Eq. (2.6) over x up to e, at fixed x + and x, we 
obtain 

-d + A-{e)+d l A*(e)- ig [ dx~ (A 1 ■ T) d + A i (x~) = f dx'UV^ p B (x~). 

Jo Jo 

(2.34) 

Replacing in this equation the field A % by its explicit expression (2.33) and per- 
forming simple manipulations (using in particular Eq. (2.16)), one can rewrite 
this equation as 

d + A-(e) = {VU)A l B (e) 

+U (W B (e) - fdx-(A B ■ T) d + A B (x-) - jT dx~ V^p B (aT)) . 

(2.35) 

The terms between parentheses that multiply U in the above equation, vanish 
since A % B solves the Yang-Mills equation (2.6) in the absence of the nucleus A. 
Therefore, we are left with 

d+A-{x) = (d l U(x+,x)) A B ( X -,x), x~<e. (2.36) 

We show now that the exact solutions (2.33) and (2.36) that are valid in 
the strip < x~ < e, extend to the region < x + < e. To do that, we 
consider Eq. (2.8) again, but now for x + < e, and recall that, in the limit e — > 0, 
$ A (x + ,x) ~ l/e. By keeping the dominant terms in the limit of small e, one 
can then reduce Eq. (2.8) to the following equation 

d-d+A l (x) -ig[<S> A (x + ,x)-T]d+A\x) = 0, (2.37) 

valid for all x + < e. This equation is easily solved and yields 

d+A*(x) = U(x+,x)d+A B (x-,x), (2.38) 

where we have used the fact that A l (x + = 0,x~,x) = A l B (x~,x). A last 
integration over x~, using A l (x~ = 0) = A l B (x~ = 0) = (see Eqs. (2.33) and 
(2.21)), gives finally 

A\x) = U(x + ,x)A B (x~,x), x + <e, (2.39) 

which is identical to (2.33), but now valid also in the strip x + < e. 

Analogous manipulations allow us to transform Eq. (2.7) into the following 
equation for d + A~(x), valid for x + < e: 

(d- - ig $ A • T) d+A- = -ig {A 1 ■ T) <9 4 $ A , (2.40) 

where we have used Eq. (2.26), (2.11) and (2.10) in order to eliminate the current 
J~. At this point we use the identity 

cT [U\x + ,x) d+A-{x)] = U*<(x + ,x) [d~ - ig $ A (x + ,x)] d + A~{x), (2.41) 
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and rewrite Eq. (2.40) as 

cT [U^x+x) d + A-{x)} = -igU\x + lX ) (A i (x) ■ T) d i $ A (x+,x). (2.42) 
By integrating over x + < e, using A~(x + = 0) = 0, one gets 

,x+ 

d + A-(x) = -igU{x + ,x) dz + U\z + ,x) {A l (z + ) ■ T) d' l <S> A {z + ,x) 

Jo 

r x+ 

= igU(x + ,x) dz + U t (z+,x)d l {<I> A (z + ) ■T)A l {z + ) 
Jo 

= ^U{x + ))A^ B {x-), (2.43) 

where in the last step we have used the expression (2.39) of A l (x), Eq. (2.43) 
is identical to Eq. (2.36), but now valid in the strip < x + < e. 

Summarizing, we have shown that the field components in the strips < 
x + < e or < x ~ < e, are given by 

A l {x) = U(x + ,x)A l B {x-,x), A-(x) = <S> A (x + ,x), A + =0, (2.44) 

and 

d+A-(x) = (d t U(x + ,x)) A l B (x-,x), (2.45) 

where A l B (x^,x) is given by Eq. (2.19). Since V(x~ > e, x) = V(x~ = e,x) = 
V(x) , A B (x~ > e,x) = A B (x) is constant in the forward light-cone x^ > e. 
The same holds for U : U (x + > e, x) = U(x + = e.x) = U {x). It follows that at 
the border of the small strip of width e along the light-cone, the field is constant 
(independent of x + and x~ ) and given by 

A\x) = U{x)A B (x), A+ = A- = , 

d+A-(x) = {d l U(x)) A B (x). (2.46) 

This solution is related to the choice of the longitudinal field (2.9) for the 
nucleus A. We have seen that by performing the gauge rotaton (2.13), which 
does not affect the light-cone gauge condition A + = 0, one obtains a transverse 
field for the nucleus A, which is analogous to the field of the nucleus B (cf. 
Eq. (2.15) and Eq. (2.21), respectively). Since, in this case A + = A~ = 0, 
the field before the collision is then similar to that in the Fock-Schwinger gauge 
x + A~ + x~ A + = 0. To make the comparison more explicit, we note that the 
field produced on the light-cone by the initial transverse field of Eqs. (2.15) 
and (2.21), is obtained from (2.44) by performing the gauge rotation (2.13) that 
eliminates A~ — $ 4 . This yields (at the border of the strip where the field 
becomes independent of x + and x~): 

A\x) = A B (x) + A\ (x) , A+ = A- = 0, 

d + A-(x) = (A i A (x)-T)A i B (x). (2.47) 

These field components on the light-cone may be compared to the corresponding 
ones obtained in the Fock-Schwinger gauge (see Eqs. (22) and (23) of Ref. [15]). 
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The transverse components, first line of Eq. (2.47), are identical to those in the 
Fock-Schwinger gauge, while the derivative d + A~ is related to the transverse 
field by a similar equation, second line of Eq. (2.47), to within a factor 1/2 (see 
also Eq. (C.120) in Appendix C where a more detailed analysis is presented). 

In the next subsections we shall consider the field in the forward light-cone, 
and present a method for solving the Yang Mills equations. We shall see that 
having transverse, pure gauge fields, for the nuclei A and B, whose supports 
extend respectively to the half planes x + > and x~ > 0, leads to spurious 
"final state interactions" in the forward light-cone, all the way till t = +oo. The 
presence of these pure gauge components greatly complicates the calculation of 
observables, and it is desirable to get rid of them. We shall see that this is 
possible by a slight redefinition of the gauge condition, which will be discussed 
in the subsection 2.5. 

2.4 The gauge field in the forward light-cone 

While it has been possible to calculate exactly the gauge field on the light-cone, 
the corrcsponging calculation in the forward light-cone remains a hard problem, 
because of the non linearity of the Yang-Mills equations. Our strategy will then 
be to construct the solution through an iterative procedure. This produces a 
solution in the form of an expansion: 

oo 

where corresponds to the non-interacting part of the gauge field, i.e, 

4»= A b> A Io)=**> 4)=°- (2-49) 
The first correction A^ is obtained by solving the linearized Yang-Mills equa- 
tions (in A^, keeping A^ as a background field) in the forward light-cone, 
using as initial conditions the field on the light-cone, which is known exactly. 
The process is repeated for A^ 2 y and so on. Note that the boundary condition 
on the border of the strip {x + = e, x~ > e; x~ = e, x + > e} reads = for 
n > 1. That is, near the light cone, the field is entirely given by A^ + A*ty 
There is therefore some analogy between the expansion (2.48) and the proper 
time expansion introduced in Ref. [29], but the two expansions arc distinct. 

In this subsection we shall construct A^ explicitly, and show that its cal- 
culation can be greatly simplified by a convenient choice of gauge. We shall 
introduce the following notation for the fields on the light-cone (Eqs. (2.46)): 

al(x) = A i (x)-A i B (x) = (U(x) (2.50) 

P (x) = d+A-(x) = (tfUixj) A* B (x). (2.51) 

Let us consider then the Yang-Mills equations in the forward light-cone, i.e., 
for x > e. The linearized equations read 

D (0),F { 7) + %)^fo) =0, (2.52) 
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where = d M — igA'Ls ■ T, and is given in Eq. (2.49) above. Since & A 
is confined to the strip < x + < e, only A l B contributes to and therefore 
d~ and D l , Q ~. = d % — igA B ■ T. Furthermore, since A % B is a 
0, and the second term vanishes. We are left with 



pure gauge, 



Making the components explicit, we get 



d+A^-D^A^ 



= 0, 



d + d-A- {1) Dl 0) Dl 0) A w + d-D^A^ = 0, 



'A: 



(2.53) 

(2.54) 
(2.55) 
0. (2.56) 

(2.57) 

where we have used the fact that d + A~(x~ = e) = /3 and A l (x~ = e) = a*. 
By using this relation (2.57) in Eq. (2.56), we eliminate the non-diagonal terms, 
i.e., the last two terms in the left hand side; we get then, after some algebra 



2d + d A l {1) D^D^A^ + D l (0) D J {0) A J {1) - D l (a) d^ 
By integrating Eq. (2.54) over x~ , from e to x~ > e we obtain 



(2d+d~ - Df 0) ) A\ x) = D\ 0) (fi D> oft 



(2.58) 



where D 



(0) 



d 1 r> % 



In order to solve this equation, we introduce the Green's function G B (x, y) 
that describes the propagation of a gluon in the background field A l B , that is, 
G B (x,y) is solution of the equation 



(2d+8- - D 2 {0) ) G B (x,y) = S(x - y) 



(2.59) 



with retarded conditions (G B is non- vanishing only when x ± > y^). The gauge 
field can then be expressed as follows (see Appendix B for details of the deriva- 
tion) 



A^x) 



d 2 y G B (x,y) a*(y) 



d A y G B {x,y) 



D 



(0) 



(3 {y) - D\ Q) <iy) 



(2.60) 



y+,y >e 



It remains to solve Eq. (2.59) for G B (x,y). To this aim, we note that for 
any color vector F, we have 



D\ 0) F = V^d\VF) and D 2 (0) F = V ] d 2 {VF) 



(2.61) 
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Thus, with □ = 2d + d — D%\, we can write Eq. (2.59) as 

V^n(VG B (x,y))=S(x-y), (2.62) 

or [28] 

G B (x,y) = V^(x)G (x-y)V(y). (2.63) 

We used the fact that V is independent of x + and x~ in the forward light-cone, 
and G (x — y), the free retarded propagator is known [24]: 

G (x -y) = 9{x+ - y + )9(x- - y~)8{{x - y) 2 ). (2.64) 

The Wilson line V attached to the two ends of the propagator G B in Eq. (2.63) 
reflects the presence of the pure gauge field A l B in the forward light cone. Would 
it be absent, Gb would reduce simply to the free retarded propagator. Indeed 
we shall see shortly that a simple redefinition of the gauge condition allows us 
to eliminate A g , thereby simplifying the calculation. 

Before doing so, it is useful to transform Eq. (2.60) into a simpler equation. 
This is obtained by applying the operator DV on both sides of Eq. (2.60). Then, 
by using the fact that D(VG B ) = V, one gets 

U{VA l (1) ) = 26(x + )5(x-)Val + 6(x+)6(x-)d l (V ' p a - d J {Va^)), (2.65) 

where we have taken the limit e — > 0. We note that the fields ^m, a* and (3 all 
appear multiplied by the same Wilson line V, which can therefore be absorbed 
in a redefinition of the fields: A 1 ^-. — VA 1 ,^, a l o — Va % and (3 = Vf3 . The 

equation for the gauge field A^ is then simply: 

□Afo = 2S(x+)S{x-)ai +8{x+)8{x-)d l o - (2.66) 

This redefinition of the fields involves the same gauge transformation as that 
which relates Gb and Go in Eq. (2.63). This suggests that the calculation 
can indeed be made simpler by a proper redefinition of the light-cone gauge 
condition. This is discussed in the next sub-section. 

2.5 The light-cone gauge d~A + = 

Consider then the gauge rotation induced by the Wilson line V: 

A> ■ T = VtA* ■ T)V ] Vd^vK (2.67) 

The effect of this gauge rotation on the field A^ is given by 

A+yT = -1-Vd + V^<S> B -T, 

A7 oy T = V{A W -T)V^ --Vd-V^V{^ A -T)V\ 

A\ 0) ■ T = V(A\ Q) ■ T)0 - — Vd l V f = 0. (2.68) 
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This gauge rotation eliminates the transverse component A l B and generates a 
non-zero A + component. Strictly speaking this gauge held is therefore no longer 
compatible with the light-cone gauge A + — 0. However, nothing prevents us 
to redefine the gauge condition so that A + = <5 B . Since <1> B exists only in a 
small strip < x~ < e, this is equivalent to the condition A + = almost 
everywhere, in particular in the forward light-cone. In fact, the slightly more 
general definition of light-cone gauge, 

d~A + = 0, (2.69) 

includes all the gauge choices that we have considered previously, including of 
course A + = and A + = <I> B , and enlarges the possibilities for fixing the initial 
fields. 

Consider now the effect of the gauge rotation (2.67) on the higher order terms 
of the expansion (2.48). The inhomogencous term in Eq. (2.67) is already taken 
into account in the transformation of into A^, Eq. (2.68). Therefore, for 
all n > 0, the fields A^ transform homogeneously 

A{ n yT = V(A{ n yT)Vl, (2.70) 
or, by using the identity VT a V^ = T b V ba , 

if , = VAf , . (2.71) 

(n) (n) V I 

This applies in particular to A^, and this completes the proof that the redef- 
inition of the fields from Eq. (2.65) to Eq. (2.66) is indeed a gauge rotation 
involving the Wilson line V. 

After performing the gauge rotation (2.68), which removes the pure gauge, 
the zeroth order field has the following components 

A to) = **> A Io)= V ®a, 4>)=0, ( 2 - 72 ) 

with the non vanishing components entirely localized in the small strip of width 
e along the light cone (see Fig. 3). The gauge field at the border of this strip 
are then simply those of Eqs. (2.50) and (2.51), to within the gauge rotation 
induced by V: 

&i=Vai=V(U-l)A i B , (2.73) 
p =Vp = (&U)Ai, (2.74) 

and are produced by the following conserved current (see Eq. (2.30)) 

J" = 6»-Vp A + S^+VUV^Pz. (2.75) 

The particular choice of gauge discussed in this subsection makes the calcu- 
lation of gluon production easier and physically more transparent since all the 
pure gauge fields are eliminated in the forward light-cone . 
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Nucleus B 



Nucleus A 



Figure 3: The gauge field before the collision in the light-cone gauge A + = <5 B 
is localized in a small strip of width e along the light-cone. 

3 Gluon production 

We now exploit the simplification brought by the use of the modified gauge 
condition A + = $ B , in calculating the spectrum of gluons produced just after 
the collision. As argued earlier only the leading correction in the expansion 
(2.48) is needed in this calculation, that is, -^m- Let us recall that the zeroth 
order A^ is given by Eq. (2.72). Since in the forward light cone x^ > e, 
$ 4 = 5> B =0, we have A^ = 0. It follows that in the forward light cone, 
D) Q \ = d l , and the linearized Yang-Mills equations take the simple form: 



Using the same technique as in the previous section to handle the boundary 
conditions, we get from them the following equations for the non vanishing 
components of the gauge field: 



As expected, the first equation (3.79) is identical to Eq. (2.66). As we shall see, 
this is in fact the only one needed to calculate gluon production. 



d + (d^ {1) )=0, 



(3.76) 
(3.77) 
(3.78) 



= 26{x+)S(x-)al + 6{x+)O{x-)d l o - d j &i), 
□ir i} = 25{x + )6{x-)f3 a . 



(3.79) 
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The spectrum of produced gluon is given by (see e.g. [24]) 

^f2E^- = Y J {\Mx\ 2 ), (3-80) 

where q is the three-momentum of the gluon. E = |q| its energy, A its polariza- 
tion. The symbol (...) stands for the average over the color sources p A and p B . 
The amplitude ftA\ is related to the classical gauge held by 

Mx = hm q 2 A l (q)e\(q), (3.81) 

where e\ (q) is the polarization vector of the gluon and q its four- momentum. In 
the axial gauge A + = 0, only the transverse components of the held contribute, 
and the sum over polarizations states is done with the help of the relation 

£4(<7K J '(<7) = ^'- 

A 

It is easy to see that only A l ^(q) contribute to gluon production (A l ^{q) 
vanishes in the forward light-cone). In Fourier space Eqs. (3.79) read 

~1 2 ^) = -2 - Sib) - «2^A(9), 

-q 2 A^(q) = (3-82) 

with 

6P Q {q) = J (fxe-^&Kx), (3.83) 

and similarly for (3o(q). Note that, for on-shell gluons (2q + q~ — q 2 = 0), the con- 
dition of transversality q il M.^ 1 {q) = is fulfilled, where ./VP = lim 9 2^ q 2 A^(q). 
By inserting in (3.81) the explicit expression of A^(q) given in Eq. (3.82) we 
get 

M X = -2 (4 - &i(q) 2i^8 ( Q ). (3.84) 

At this point we make the following decomposition of a* 

ai^^(q-a (q)) + e l ^(qxa (q)), (3.85) 

where e 1 -* is the antisymmetric tensor with e 12 = 1. This allows us to write the 
gluon spectrum in the following compact form: 

47r3 ^ = Y ilq X + ^° (9)|2> ' (3 ' 86) 

The formula (3.86) is the main result of this paper. It provides a compact 
expression for the spectrum of gluons produced in a nucleus-nucleus collision 
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in terms of the fields a and $ that exist on the light-cone just after the 
collision. As shown by their explicit expressions (2.73) and (2.74), these fields 
are determined entirely by the functions <1> A and <f> B that are obtained as a 
function of the sources p A and p B by solving Eqs. (2.10) and (2.18). The complex 
color structure of the product of the Wilson lines U and V hidden in a and /3 
prevents one to write Eq. (3.86) in a fc f -factorized form. This property of k t - 
factorization is recovered, however, when one of the sources is weak, as we shall 
see shortly. The spectrum (3.86) is rederived in Appendix C using a variant of 
the Fock-Schwinger gauge which makes the calculation more symmetrical: this 
provides an indication of the robustness of the calculation against variations in 
the gauge choice. 

The formula (3.86) is obviously not exact since it is the leading order of 
an expansion (it involves only -^n))- However the initial values of the fields 
on the light-cone are calculated exactly, i.e., taking the full non-linearity of 
the Yang-Mills equations into account. Such non-linearities are ignored in the 
calculation of the field Am which propagates in the forward light-cone without 
further interactions. The next order in the expansion (2.48), A( 2 ), describes in 
particular the merging of such gluons. Note that such mergings are taking place 
at later times since A( 2 ) = when x ± = e. 

The formula (3.86) becomes exact when one of the nuclei can be considered 
as a weak source. To see that let us expand the Wilson U line at leading order 
in P A 

U(x) = 1 + ig [ dx + $ A (x + ,x) • T, (3.87) 
Jo 

which allows us to write 

al(x) = igVT a A^ a A {x) = {&V)* A {x), (3.88) 

and 

o {x) = igVT a A j B d^ a A {x) = {d j V)d' J <$> A {x), (3.89) 
where we have used the identity 

(VT a A l B ) b = -{d*V) bai (3.90) 
ig 

and we have set dx + $ A {x + 1 x) = <fr A [x). In Fourier space, setting p = q — k, 
we have 

= -if V( P ) <Mfc), (3.91) 

0M = ~j -^( P -k)V(p)$ A (k), (3.92) 



with 



(27T) 2 

V(p) = J d 2 x e- ip - x (V(x) - 1). (3.93) 
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By replacing in (3.86) 6t l Q (q) and (3 (q) by their expressions (3.91) and (3.92) one 
recovers the well known formula for proton-nucleus collisions (see e.g. Eq. (22) 
of Rcf. [25]): 

^ E ik = ?/ w?w? /c(fe)fe ' q) {v ' {p ' )v{p)U wxwk 

(3.94) 

where p' = q — k' and JC(k,k' ,q) is the square of the Lipatov vertex 

K(k,k',q) = (px fe)(p' x fc') + (p-fe)(p'-fc')- (3-95) 

This formula exhibits fc t — factorization. 

Although the formalism that we are using is not manifestly symmetric (be- 
cause of the gauge choice), it is easy to verify that a similar result can be 
obtained in the case where the weak source is p B : 

K(q) = -iJ -^P 1 U(k) $ B (p), (3.96) 

P°( q) = ~ j ' -^(P- k )U(k)<P B (p), (3.97) 

which leads to 

(3.98) 

This is identical to Eq. (3.94) after the substitutions A — > B, V — > U and 
p <-> k, p' <-> fc', under which the kernel /C is invariant. 

4 Conclusion and perspectives 

We have shown that the calculation of gluon production in nucleus-nucleus col- 
lisions in the color glass framework can be greatly simplified by an appropriate 
choice of gauge. The particular axial gauge that we have proposed, d~A + = 
(or more precisely A + = $ B ), allows us to determine explicitly the conserved 
current, and the fields immediately after the collisions, in terms of the color 
sources of the two nuclei, given in the covariant gauge. The fields before the 
collision are localized in an infinitesimal strip along the light cone. In particular, 
no pure gauge fields are left in the forward light cone; in other gauges the pres- 
ence of such pure gauges complicates the calculations and sometimes obscures 
the physical interpretation. 

Wc have obtained an explicit expression, Eq. (3.86), for the spectrum of glu- 
ons produced in nucleus-nucleus collisions, which holds in leading order in a new 
expansion where the non linearities of the Yang-Mills equations in the forward 
light-cone are treated itcrativcly. The leading order is obtained by linearizing 
the Yang-Mills equations in the forward light cone, keeping as boundary values 
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for the fields, the exact fields that exist on the light-cone immediately after the 
collision. Higher order terms in the expansion describe "final state interactions" , 
in particular the merging of gluons; these, however, take place at later times, so 
that our formula exactly calculates the spectrum of gluons produced at a short 
time after the collision. It would be interesting to investigate quantitatively 
the importance of the final state interactions, for instance by comparing (3.86) 
to available numerical studies for gluon production [13,14] and for the energy 
density [13,30], and check whether these are negligible, as argued in Ref. [18]. 

The iterative solution that we propose is not an expansion in an obvious 
small parameter. However it lends itself to simple calculations whenever such a 
parameter is present. Thus, for instance, the formula (3.86) provides the exact 
leading order result when expanded in powers of the color source of one of the 
two nuclei. An immediate application of the present formalism is to obtain the 
ncxt-to-lcading order in this expansion. This requires the determination of the 
term A 1 ^ in our expansion. The result of this calculation will be presented in 
a forthcoming publication. 



A A along the x + axis 

In this appendix we establish the continuity of A~(x~) at x~ = 0, that is, we 
show that 

\mxA~(x + ,x~ = e,x) = $ A (x + ,x). (A. 99) 

e— »0 

To do so, let us rewrite the Yang-Mills equation (2.6) in the form 

d + (d + A~) = d + {d l A l ) - J+ - ig{A l ■ T) d + A\ (A.100) 

and integrate it over x~ from — oo to x~ . Using the initial conditions A m (.t _ < 
Q 1 x + ,x) = d>- L ~ $ A (x + ,x) and J + (x~ < 0) = 0, we obtain 

d+A-( x -) = 0^4* (aT) " / dy~ [J + (V~) + ig(A< ■ T) d+A l (y-)] , (A.101) 

Jo 

where the dependence on x + and x has been omitted to alleviate the notations. 
After a further integration of Eq. (A. 101) over x~ , from to e, one gets 

A-(e)-$ A = f dx-d l A l (x~) 
Jo 

- f dx~ f dy- [j + {y-) + ig{A l ■ T) d+A l (y-)]. (A.102) 
Jo Jo 

The continuity of A~{x~) follows then from the absence, in the limit e — > 0, 
of 5-function singularities in the integrands in the right hand side (such singu- 
larities occur in the current J + , or in d + A l , but these disappear after the first 
integration, over y~). 
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B Proof of equation (2.60) 

We are looking for a solution of an equation of the form 

(2d + d--Df 0) )A i (y) = J i (y), (B.103) 

in the region x + > e and x~~ > e, which we shall obtain by using the retarded 
Green's function G B defined in Eq. (2.59). Let us multiply Eq. (B.103) by G B 
and integrate over y in the forward light cone: 

d 4 yG B (x iy )(2d + d-~ Df 0) )A i (y)= J d±y G B (x lV )J\y), 

y + ,y~>t y + ,y~>£ 

(B.104) 

Loosely speaking, the Green's function G B is the inverse of the differential 
operator which sits next to it in Eq. (B.104). However, care must be exerted in 
using this property because the y-integration covers only the forward light-cone 
and, as we shall see, boundary terms will survive integrations by parts. 

Let us first consider the first term in the left hand side of Eq. (B.104). 
Performing successive partial integrations we end up with 



d A y G B (x, y) 8+d- A 1 (y) = J d'y G B (x, y) d+d~ A 1 (y) 
J dy + d 2 yG B {x,y)d~ A l (y)+ J dy~ d 2 y G B (x,y) d+ A <(y) 

y+=e 

d 2 yG B (x,y) A l (y), (B.105) 



y T ,y =e 



Since A l (y + — e, y~ > e,y) = A l (y + > e,y~ = e,y) = A l (y), as shown in Sect. 
2.3, the integrals over y + and y~ in the second and third term in the right hand 
side can be performed explicitly, leading to 



dS G B (x, y) dp- A 1 (y) = J d+y G B (x, y) 8+8; A { (y) 
y + ,y~>£ y + ,y~>£ 

- J d 2 yG B {x,y)A\y). (B.106) 

We turn now to the second term in the left hand side of (B.104). By performing 
successive partial integrations, and assuming that the field A l {x) vanishes when 
\x\ — ► oo we get 

J d A y G B (x,y) D 2 {0)y A\y) = J d 4 y G B (x,y) E>\ 0)y A l (y), 

(B.107) 
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y + ,y >e 



where DU = d l + igA l B ■ T. Finally, combining Eqs. (B.106) and (B.107), we 
obtain 

r <- <- <- 2 

A G B (x,y)(2 fl+a- - Dt (0)5/ ) A%) = A\x) 

= 2 / d 2 yG B (x,2/)^(y)+ / d±y G B (x,y) J\y) 7 

y + ,y~=£ y + ,y~>e 

(B.108) 

from which Eq. (2.60) follows, with J 1 = D| 0) (/3 - D 3 (0) a^. 

C Comparison with the Fock-Schwinger gauge 

In this appendix we briefly review the calculation of the gauge field in the axial 
gauge 

x + A~ + x~A + = 0, (C.109) 

often referred to as the Fock-Schwinger gauge, and compare with the results 
obtained in this paper. In the gauge (C.109), we have A~ = along the x + 
axis, and A + = along the x~ axis. It is then easy to verify that the current 

j+ = v^ Pb , j- = u^ Pa , r = o, (c.iio) 

is covariantly conserved, and that the gauge field components before the collision 
are transverse: 

A 1 = A\ for x+ > 0, 
A 1 = A l B for x~ > 0, 

a + = a~ = o, (cm) 

where 

A\-T= -l-Ui&U, A\-T= -—V^d'V, (C.112) 
ig ig 

with U = U(x+,x) given by Eq. (2.14) and V = V(x~,x) by Eq. (2.20). 

The gauge choice (C.109) suggests the following ansatz for the components 
A ± (x) of the gauge field after the collisions [15]: 

A + {x) = -x + P fs (t,x), A~{x) =x-P fs {t,x), (C.113) 



where r = 2V x + x is the proper time. The gauge field A(x) in the light-cone 
gauge is related to the gauge field A{x) of the Fock-Schwinger gauge by a gauge 
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rotation: 



a+-t = n\A+ -T)n- —rtd+n = 0, 



A" 



T = n^(A~ -T)Q- — Q^d-fl, 



A 1 ■ T = at(A i -T)Sl-—ti i d i {l. 
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(C.114) 



Is is easy to verify that, given the ansatz (C.113), there exists a boost invariant 
solution of the equations above, that is, a function £1 that depends only on 
the proper time r, £l(x + ,x~,x) = f2(r, x). To show this, we divide the first 
equation by x + and use Eq. (C.113) in order to express A + in terms of (3 FS . 
Assuming then that Q depends only on r, we obtain 



4 Afi( T) x ) + lg (f3 FS ( T , x) ■ T) n( T> x) = 0, 



whose solution reads 

f2(r, a;) = Pexp 



4 



dea FS (Lx)-T 



(C.115) 



(C.116) 



By dividing the second equation (C.114) by x and assuming for A (x) a form 
analogous to that of A~ (x) in Eq. (C.113), i.e., A~ (x) = x~j3 LO (r, x), we obtain 

(3 LC -T = tf{fi PS -T)Q-ifit_L fi 

ig ot z 

1 



nt(/j FS . t) n - -nt(_ i5/?ps . T ) n 

2Qt(/3 PS .T)fi. 



(C.117) 



Near the light-cone, i.e., when r — > 0, Q — * 1. Then, from (C.114) and (C.117) 
we get, respectively, 

A i (0,x) = A i (0,x), p FS (0,x) = ±f3 LG (0,x), (C.118) 

where we have assumed that the transverse components of the gauge field, A 1 
and A 1 , depend only on r (with *4(0, x) = A(t = 0, x) and similarly for A(0, x)). 

In the light-cone gauge, with the further choice that leads to the same initial 
(transverse) fields as in Eq. (C.112), i.e., A\ = A\ and A l B — A l B , the gauge 
field on the light-cone is given by Eq. (2.47): 



d+A- = (A* ■ T)A\ 



A 1 = A\ + A\ 



(C.119) 



Therefore, we reproduce the well known result for the gauge field on the light- 
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cone in the Fock-Schwinger gauge, namely [15-17] 1 : 
A i (0,x) = A i A (x) + A i B (x), 

d+A-(0,x) = -d-A + (0,x)=(3 F3 (0,x) = ^{A\(x)-T)A B {x). 

(C.120) 

At this point we depart from the traditional treatments in the Fock-Schwinger 
gauge. In complete analogy with what we did in Sect. 2, we eliminate the trans- 
verse pure gauge fields by performing a gauge rotation involving the products 
UV or VU of the Wilson lines (2.14) and (2.20) 2 . By doing so, we generate a 
longitudinal initial field (for t < 0) 

A + = <$> B , A~=$ A , -4* = 0. (C.121) 

Such a field is no longer compatible with the gauge fixing condition (C.109). 
Consider however the following extension of this condition, 

d~A + +d + A-=0, (C.122) 

which may be viewed as the symmetric version of the gauge condition (2.69). 
The initial field (C.121) clearly fulfills this condition. Besides, in the forward 
light-cone, and provided the ansatz (C.113) is verified, the gauge conditions 
(C.122) and (C.109) are simultaneously verified. 

We shall use now this symmetric gauge condition to calculate gluon produc- 
tion in leading order in the expansion (2.48), and recover the results of Sect. 3. 
Following the steps described in Sect. 3, one obtains first the linearized Yang- 
Mills equations in the forward light-cone (compare with Eqs. (3.76), (3.77), 
(3.78)): 

o^ij = -d-&A\ iy 

UA\ X) = -~d l d J A J {1) , (C.123) 

with A(i)(t = 0,x) = A(t = 0,x) and A{ n ){r = 0,x) = for n > 1. The 
values of the fields, created by the initial fields (C.121) near the light-cone, are 
related to those in light-cone gauge (given in Eqs. (2.73) and (2.74)), through 
Eq. (C.118). More explicitly, A 1 {t = 0,x) = a* (a:) and d + A~(T = 0,x) = 
—d~A + (r = 0, x) = ^/3 (x). The first two equations in (C.123), together with 
the gauge condition (C.122), lead to 

= 2d + d~d j A j (1) - 0, (C.124) 

the references quoted here the fields are expressed in matrix representation, i.e. A-T — > 
A and (3 FS = \[A\,A B ] 

2 The elimination of the pure gauge components of the gauge fields is useful in view of the 
calculation of gluon production, as was also observed in [23]. 



□ 



d+A 



(i) 



d~A 



(i) 
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from which it follows that the divergence of the transverse field is conserved, i.e. 

d 3 A 3 {1) =d 3 al (C.125) 

This allows us to rewrite the Yang-Mills equations in a form that makes the 
boundary conditions explicit. Following the same steps as in Sect. 2.4, we 
obtain 

04 ? (1) = 25{x+)S(x~)a l o - 0{x + )0{ X -)d l (d 3 a 3 ) : (C.126) 
04+) = S{x-)6(x + )p , (C.127) 
04 a) = 5(x + )6(x-)(3 . (C.128) 

In Fourier space we get 

-q 2 A^(q) = 2 (> - 0-} &i(q) - ±fa q ) 6* ±0 o {q) 6^. 

(C.129) 

From here, it is straightforward to get the gluon spectrum using the decompo- 
sition (3.85): 

= ±{\q x a (q)\ 2 + \P (q)\ 2 ), (C.130) 
which is identical to Eq. (3.86) obtained in light-cone gauge. 
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